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ABSTRACT 


A  demonstration  is  given  that  for  two-electron  systems  the  Hartree- 
Fock  and  first  natural  spin  orbitals  are  different  functions.  The 
method  used  is  perturbation  theory  for  which  the  zero-order  problem  is 
the  Hartree-Fock  approximation.  A  perturbation  expansion  through  the 
second  order  is  obtained  for  the  first  natural  spin  orbital.  The  two 
orbitals  begin  to  differ  in  the  second  order  and  their  energies  in  the 
fourth  order.  An  equation  for  the  second  order  part  of  the  orbital 
difference  function  is  derived.  Estimates  of  the  norms  of  the  orbital 
difference  functions  are  calculated  for  the  ground  states  of  the  helium 
atom  and  the  hydrogen  molecule  and  are  found  to  be  small.  . 
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•r'  COMPARISON  OF  THE  HARTREE-FOCK  ORBITAL  WITH  THE  FIRST 

NATURAL  SPIN  ORBITAL  FOR  TWO-ELECTRON  SYSTEMS 

Introduction 

Calculations  of  approximate  energies  of  two-electron  systems  have 
shown  that  the  first  natural  spin  orbital  gives  an  approximate  energy 
very  close  to  the  Hartree-Fock  energ^f^^-*^'^^  The  question  has  arisen 
whether  the  first  natural  spin  orbital  would  coincide  with  the  Hartree- 
Fock  orbital  If  It  were  calculated  exactly.  One  would  expect  the  two 
orbitals  to  be  different  functions  because  they  are  defined  differently. 
The  Hartree-Fock  orbital  Is  defined  to  give  the  optimum  energy^  whereas 
the  first  natural  spin  orbital  Is  defined  to  have  the  optimum  overlap 
with  the  true  wave  function.  In  this  paper  a  demonstration  Is  given 
that  the  first  natural  spin  orbital  Is  different  from  the  Hartree-Fock 
orbital  for  two-electron  systems  even  though  the  difference  is  small. 

The  method  followed  in  this  paper  is  perturbation  theory  in  which 
the  zero-order  problem  is  the  Hartree-Fock  approximation.  The  discus¬ 
sion  starts  with  a  brief  summary  of  the  proof.  After  that  follows  a 
short  review  of  natural  spin  orbitals  for  two-electron  systems.  The 
perturbation  scheme  Is  then  defined,  and  an  expression  for  the  first 
natural  spin  orbital  Is  developed  in  terms  of  the  Hartree-Fock  orbital. 
The  difference  function  between  the  two  orbitals  is  then  discussed,  and 
Its  norm  Is  explicitly  calculated  for  the  helium  atom. 

Throughout  the  discussion  the  exact  wave  function  of  the  two- 
electron  system  Is  always  the  spatial  part  of  the  singlet  ground  state 
wave  function.  Likewise,  only  the  spatial  parts  of  the  Hartree-Fock 
and  natural  spin  orbitals  are  considered.  The  spin  parts  are  disre¬ 
garded  since  the  Hamiltonian  does  not  contain  any  spin  operators. 


H.  Shull  and  P.-O.  LOwdin,  J.  C.  P.  23,  1565  (1955). 

P.-O,  LOwdln  and  H.  Shull,  Phys.  Rev.  m,  1730  (1956). 

(3) 

'■  ’  H.  Shull  and  P.-O.  LBwdin,  J.  C.  P.  30,  617  (1959). 

(4) 

'  E.  R.  Davidson,  J.  C.  P.  37,  2966  (1962). 
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In  brief,  the  proof  that  the  first  natural  spin  orbital  is  differ¬ 
ent  from  the  Hartree-Fock  orbital  is  as  follows.  The  exact  wave  function 
of  the  two-electron  system  is  first  expanded  in  a  perturbation  series, 
eq.  (12),  in  which  the  leading  term  is  the  Hartree-Fock  function.  Then 
the  perturbation  expansion  of  the  first  natural  spin  orbital,  eq.  (14), 
correct  through  the  second  order  in  the  perturbation  parameter^  is 
derived.  It  becomes  apparent,  eq.  (18),  that  the  difference  between  the 
two  orbitals  is  not  everjrwhere  zero.  Therefore  the  two  orbitals  are 
different  functions. 

The  use  of  perturbation  theory  is  justified  by  the  smallness  of  the 
perturbation.  An  order  of  magnitude  estimate  of  the  perturbation  can  be 
obtained  from  the  following  considerations.  The  difference  between  the 
Hartree-Fock  energy  of  a  two-electron  system  and  the  approximate  energy 
given  by  the  first  natural  spin  orbital  should  be  of  the  fourth  order 

in  the  perturbation^^  Shull  and  LBwdin^^^  found  this  energy  difference 

-4  2  (4) 

to  be  1.43  X  10  e  /a  for  the  helium  atom;  Davidson  obtained 

-4  2  ° 

1.57  X  10  e  /a^  for  the  hydrogen  molecule.  Therefore,  the  pertur¬ 
bation  is  small  enough  to  justify  a  perturbation  expansion. 


Perturbation  Expansion  of  the  First  Natural  Spin  Orbital 

In  this  section,  a  few  properties  of  the  natural  spin  orbitals  are 
first  discussed.  Then  the  perturbation  scheme  which  is  used  in  develop¬ 
ing  the  exact  wave  function  of  a  two-electron  system  is  explicitly 
stated.  After  that,  the  perturbation  expansion  of  the  first  natural 
spin  orbital  is  developed. 

The  exact  normalized  wave  function  i|r(12)  of  a  two-electron  system 
can  be  expanded  in  terms  of  the  natural  spin  orbital^^'’^'^'^^  ^ as 
follows: 


See  footnote  (8). 

P.-O.  LBwdln,  Phys.  Rev.  1474  (1955). 
E.  R.  Davidson,  J.  C.  P.  37,  577  (1962). 
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i|r(12) 


•D 

I’l 

k-1 


(2) 


(1) 


where  the  are  the  occupation  numbers.  The  natural  spin  orbitals 

are  mutually  orthonormal,  l.e. 


«1  -  . 


(2) 


For  the  purposes  of  this  discussion  It  Is  convenient  to  rewrite  the 
wave  function  iir(12)  as 


♦  (12)  - 


+  T(12j 


(3) 


where  T(12)  contains  all  of  the  natural  spin  orbitals  except  the  first 
one.  An  Important  property  of  the  first  natural  spin  orbital  Is  that 
It  Is  strongly  orthogonal  to  the  function  T(12)  ,  l.e. 


jT*(n)%^a)d\  -  ^*(12)X^(2)  dt^  -  0  . 


(4) 


The  perturbation  scheme  Is  obtained  by  rewriting  the  exact 
Hamiltonian  H  of  the  two-electron  system  as  a  Hartree-Fock  Hamilton¬ 
ian  H^°^  plus  a  perturbation  H^^^ 


H  -  H<°>+>H<^> 


(5) 


where  \  Is  the  perturbation  parameter  which  Is  ultimately  set  equal 
to  unity.  The  Hartree-Fock  Hamiltonian  is 


i(o) 


h^d)  +  h‘’(2) 


(6) 


where 


h°(l;  -  vf  -  -r  +  V(l) 

1 


(7) 


4 

In  atomic  units.  The  Hartree-Fock  potential  V(l)  In  which  electron 
1  moves  Is 

The  function  is  the  Hartree-Fock  orbital  which  satisfies  the 

equation 


(9) 


with  as  the  orbital  energy.  The  Hartree-Fock  orbital  is  normalized 

to  unity,  i.e. 


(1)  dCj  -  1 


(10) 


The  Hartree-Fock  function  of  the  two-electron  system  is  given  by  the 
product  fc  (2)  .  The  perturbation  term  H^^^  in  the  Hamiltonian 

is 


H^^^  =  -L.  .  .  v(2) 

*^12 


(11) 


where  V  is  defined  by  eq.  (8),  and  rj^2  1®  inter-electronic 
distance. 

According  to  the  above  perturbation  scheme  the  exact  wave  function 
i|f(12)  is  written  as  yg(l)^p(2)  plus  additional  terms.  The  detailed 
form  of  this  expansion  through  the  second  order  in  A  is  derived  in 
Appendix  1,  eq.  (1,10),  and  is 


y(12)  -  c^^(l)f^(2)  +  jy^(l)f^2\2)  +  f^^\l) 

+  >S<^\l2)  +  A^<2>(12)  + 


(12) 


The  normalization  constant  C  is  also  the  overlap  between  y(12)  and 
the  Hartree-Fock  function  ^j,(l)yp(2)  ■  ^1*®  functions  f^^^  , 


obey  the  following  orthognallty  relations  (see  end  of 


and  S<2) 
Appendix  1) , 


f 


dr^ 


(13a) 


and 


/.a.*, 


/ 


(12)y^(l)dc^  -  0 


s(2)(i2)y^(l)dt^ 


(13b) 


(13c) 


From  an  examination  of  ^(12)  given  by  eq.  (12)  It  is  seen  that 
the  second  order  expansions  for  the  first  natural  spin  orbital  and  the 
function  T(12)  appearing  in  eq.  (3)  are 


and 


X^d)  -  ^^(1)  +  AV^^I) 


T(12)  -  As^^^(12)  +A^S^^\i2)  . 


(14) 


(15) 


The  overlap  between  9  (1)  and  X.(l)  is  equal  to  unity,  through 
Jo  i  ,2\ 

the  second  order,  since  ^^(1)  is  orthogonal  to  f'  '(1)  ,  eq.  (13a). 
Therefore,  the  overlap  C  between  +(12)  and  yQd)^p(2)  is  approx¬ 
imately  equal  to  ,  the  square  root  of  the- ft  -,  t  occupation  number. 

The  two  following  observations  indicate  that  eqs.  (14)  and  (15) 
are  the  correct  second  order  expansions.  First,  substituting  the  two 
expansions  into  +(12)  given  by  eq.  (3)  and  keeping  terms  only  through 
the  second  order  in  ^  gives  the  correct  form  of  +(12)  through  the 
second  order,  eq.  (12).  Second,  the  perturbation  expansion  of  the 
first  natural  spin  orbital,  eq.  (14),  satisfies  the  normalization 
condition,  eq.  (2)  with  k  »  k'  ■  1  ,  through  the  second  order  in  A  , 
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by  virtue  of  eqs.  (10)  and  (13a).  The  two  expansions  also  satisfy^ 
through  the  second  order,  the  strong  orthogonality  condition,  eq.  (4) , 
l.e. 

y'[\s^^^l2)  +)(^S^^\l2)J  jy^(l)  dt^  -  0  (17) 

by  virtue  of  eqs.  (13b)  and  (13c).  Therefore,  it  can  be  concluded  that 

eq.  (14)  is  the  correct  second  order  expansion  of  the  first  natural 

spin  orbital  for  a  two-electron  system. 

From  eq.  (14)  it  is  apparent  that  the  first  natural  spin  orbital 

(,  is  different  from  the  Hartree-Fock  orbital  .  The  difference 

between  the  two  orbitals,  up  through  the  second  order  in  ^  ^  fs  given 

(81  (2) 

by  the  difference  function'  f'  '  . 


(2) 

The  Difference  Function  f 

(2) 

The  equation  which  the  difference  function  f  satisfies  (see 
Appendix  II,  eq.  (11,6))  is 


(h'd)  - 


where  the  function 


.(1) 


satisfies  the  equation 


(18) 


(19) 


The  quantities  and  are  the  first  and  second  order  energies 

of  the  two-electron  system  and  are  given  by  the  well-known  eqs.  (11,3) 
and  (11,5). 

From  eq.  (19)  it  is  apparent  that  is  not  everywhere  zero. 

Therefore  it  follows  from  eq.  (18)  that  the  difference  function  f^^^ 
is  not  everywhere  zero.  Thus,  the  first  natural  spin  orbital  begins 
to  differ  from  the  Hartree-Fock  orbital  in  the  second  order. 


r  8^ 

'  From  eq.  (14)  one  can  verify  that  the  Hartree-Fock  energy  begins  to 
differ  from  the  first  natural  spin  orbital  energy  in  the  fourth 
order.  This  result  is  obtained  by  taking  the  expectation  value  of 
the  Hamiltonian  H.  ,  eq.  (5),  with  ')Lj(l)‘)ti  (2)  and  making  use  of 
the  form  of  the  perturbation  H'^'  ,  eq.  (11). 
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(2) 

An  estimate  for  £  can  be  obtained  from  a  sufficiently  good 
approximate  two-electron  wave  function  by  the  following  method;  Multi¬ 
plying  i|r(12)  ,  eq.  (12),  by  ^^(2)  ,  Integrating  over  the  coordinates 
of  electron  2  and  taking  the  orthogonality  relations,  eqs.  (13a),  (13b) 
and  (13c),  Into  account,  we  obtain 


where 


0  .^*(l2)JI^(l)JI^(2)di,d,j  . 


(20) 


(21) 


=(2) 


Equation  (20)  la  an  approximation  because  the  formula  gives  f''~''  plus 

(2) 

higher  order  contributions.  The  norm  of  f'  '(1)  In  general  Is  defined 
to  be 


p<2). 


f(2)(i)f(2)(i) 


(22) 


In  terms  of  the  two-electron  wave  function  of  eq.  (12), 


(12)JI^(2)dt2j  m  (12)J>^(2)dt 


(23) 


Equations  (21)  and  (23)  were  used  to  obtain  an  estimate  for  the 

(2) 

norm  of  f  for  the  case  of  the  helium  atom.  The  normalized  Hylleraas 
functloi^^^ 


♦  (12) 


(1.380)  \  1  +(0.3534)u+  (0.1282)t  -  (0.1007)s 


+  (0.03305)8  -  (0.03173)u  exp(-1.8188) 


(24) 


(9) 


This  method,  due  to  0.  Slnanoglu,  was  communicated  to  us  by  S. 
Tuan. 


(10) 

'  '  E.  A.  Hylleraas,  Z.  Physlk  54,  347  (1929). 
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W48  taken  as  an  approximation  to  the  exact  helium  wave  function.  The 
variables  s,  t  and  u  are  the  well-known  combinations  s  -  (r^  4-  r2)  , 
t  ■  -  rj)  and  u  »  r^^^  where  the  >  ''^2  *  *^12 

expressed  in  units  of  a^  with  being  the  reduced  mass  of 

the  helium  atom. 

The  difference  function  f^^^(l)  ,  eq.  (20)  has  the  form 


(25) 


where 


2«  r" 

‘  ^Jr  1  .‘*’^12  "2^12  ^<12)  5>^(2) 

l^»-0  r.-r..  ■' 


(26) 


12  '"1  ‘2' 


.(11) 


The  following  analytic  approximation'^ to  ^^(2)  was  used 
y^(2)  ■  (2.968466)  1.455799  r^)  +  (2.00)  exp(-0. 8734794  rjjj  .  (27) 


Here  r2  Is  In  units  of  a^  .  The  function  1(1)  turned  out  to  be 

(with  r.  In  units  of  a  ): 

1  o 

1(1)  -  1^(0.1196)  exp(-1.818  r^^)  -  (0.1092)  exp(-5.092  r^ 

-  (0.01041)  exp(-6.547  r^l  r"’^  +  (1.037)  exp(-1.818  rj^) 


-  (0.08935)  exp(-5.092  r^ 
+  (0.08650)  exp(-1.818  r^) 
The  coefficient  C  ,  eq.  (21),  is 


-  (0.01231)  exp(-6.547  r^ 

+  (0.1286)  exp(-1.818  rj^) 

0.9961  . 


(28) 

2 

’^l  • 


(29) 


(11) 


L.  C.  Green,  M.  M.  Mulder,  M.  N.  Lewis  and  J.  W.  Woll,  Phys.  Rev. 
93,  757  (1954). 
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The  norm  of  £ 


(2) 


,  eq.  (23),  la 


2 

dr^r^KDld) 


-)■ 


0.002 


(30) 


(2) 

Thus  the  norm  of  £  la  very  amell  compared  to  unity,  the  norm 
The  overlap  between  the  Hartree-Fock  and  flrat  natural  apln 
orbltala,  through  the  fourth  order,  la 


(l)Xi(l)drj 


0.999998  . 


(31) 


Our  value  for  the  overlap  la  higher  than  the  value  obtained  by  LBwdln 

and  Shull,  '  0.99995080  .  Slnano^lu  haa  mentioned  that  aa  one  Improves 

an  approximate  two-electron  wave  function,  the  magnitude  of  the  function 
(2) 

f  la  reduced.  Since  the  energy  obtained  by  LBwdln  and  Shull  la 
2  2 

-2.8785973  e  /a  aa  compared  with  -2.90324  e  /a  obtained  by  Hylleraas 
o  o  >2) 

the  wave  function  calculated  by  LBwdln  and  Shull  haa  a  larger  f  than 

the  Hylleraaa  function,  eq.  (24).  Thua,  the  above  overlaps  are  In  the 

correct  relatlonahlp. 

(4) 

Davldaon  haa  examined  the  natural  apln  orbital  expansion  of  the 

Kolos  and  Roothaan  wave  function  for  the  hydrogen  molecule.  He  obtained 

0.999982  for  the  overlap  between  the  Hartree-Fock  and  first  natural 

(2) 

spin  orbitals  of  H2  .  Therefore,  the  norm  of  the  £  associated 
with  the  Kolos  and  Roothaan  wave  function  is  approximately  0.006  . 

This  value  is  small  compared  to  unity,  the  norm  of  the  Hartree-Fock 
function  for  H2  . 
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APPENDIX  I 

Perturbation  Expansion  of  an  Exact  Two-Electron  Wave  Function 


In  this  appendix  a  derivation  la  given  of  the  form  which  an  exact 
two-electron  wave  function  ^(12)  must  necessarily  assume  If  It  Is 
developed  according  to  the  perturbation  scheme  described  by  eqs.  (5) 
through  (11). 

If  the  exact  two-electron  Hamiltonian  of  the  system  Is  decomposed 
according  to  eqs.  (5)  and  (6),  then  the  two-electron  SchrCdlnger 
equation  which  ^'(^2)  satisfies  can  be  written  as 

^h°(l)  +  h°(2)  ♦(12)  -  Bt(12)  (1,1) 

where  E  Is  the  exact  energy  of  the  two-electron  system. 

The  zero-order  two-electron  problem  associated  with  eq.  (1,1)  Is 

(h^l)  .  h°<J))j)^(l)J>,(J)  .  .  (1,2) 

Therefore  ^(12)  must  be  written  as  the  solution  of  the  zero-order 
problem,  j  other  functions. 

The  one-electron  problem  equivalent  to  eq.  (1,2)  Is 


h»(i)f„(i)  -  £„f,(i)  . 

The  one-electron  Hamiltonian  h°(l)  ,  eq.  (7),  Is 

h°(i)  .  -%  v^d)  - ^  f„(2)dtj 


(1,3) 


(1,4) 


with  the  Hartree-Fock  potential  V(l)  ,  eq.  (8),  written  out  explicitly. 

It  should  be  noted  that  the  operator  h°(l)  In  eqs.  (1,1)  to  (1,4) 

depends  on  the  unknown  function  . 

After  eq.  (1,4)  Is  solved  for  O  ,  one  can  use  the  operator 
o  I  ® 

h  (1)  with  C  fixed  to  define  the  following  eigenvalue  problem 


h°(l)f^(l)  - 


k  -  1,  2,  ... 


(1,5) 
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The  Cwo-electron  wave  function  ^(12)  can  be  expanded  in  terma  of 
products  of  the  form  where  i,J  -  0,  1,  2,  ...  .  In  these 

products  the  orbitals  are  any  solution  of  eq.  (1,5)  or  the  Hartree-Fock 
orbital  .  Therefore  the  perturbation  expansion  of  i|f(12)  through 
the  second  order  in  X  is 

(  *  *r 

♦(12)  -  cU  (l)f>,(2)  fi(l)fj(2) 


l>J-0 

higher  order  terms 


■] 


(1,6) 


where  C  is  the  normalization  constant.  The  summation  sign  includes 
Integration  over  the  continuum  while  the  asterisk  implies  that  the 
subscripts  i  and  J  are  never  both  equal  to  zero. 

The  perturbation  expansion  of  eq.  (1,1)  is  obtained  by  substituting 
into  eq.  (1,1)  from  eq.  (1,6)  and  E  from  the  well-known 

expansion 


E  -  2^  +  •••  .  (1,7 

Collecting  coefficients  of  X  snd  X  one  obtains,  respectively,  the 
well-known,expressions 


(1)  -  H‘^>|y,(i)f„(2)> 

'H  ■  (£l  +  tj  -  2  V 

and 

(2)  V 

'll  '  L  Ui  +  £,  - 

m,  n»0 

*  («.  *  ■  *  V 


where  the  brackets  are  a  short-hand  notation  for  the  integrals.  A 
straight  forward  calculation  shows  that  the  Integrals 
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one  of  the 

■ubacrlpts  1  or  J  is  equal  to  zero. 

Subatltuting  eqs.  (1,8)  and  (1,9)  into  eq.  (1,6)  one  obtains  the 
desired  second  order  expansion  of  \|r(12)  ,  namely 


♦  (12)  -  c|f^(l)^^(2)  + 

■f<^>(l) 

f„(2)  + 

+  Xs^^\l2)  + 

y,i2) 

(1!)  + 

The  functions  f^^^  ,  and 

s(2> 

are 

.  1  c<=).J 

1-1 

»i(l) 

> 

CO 

(l)f^^\2)! 


terms 


•1 


S^^^12)  -  ^  , 


and 


S^^^12)  -  ^  • 


l,j-l 


(I. 10) 


(I,lla) 


(I, lib) 


(I,llc) 


Since  the  Hartree-Fock  orbital  ^  Is  orthogonal  to  all  of  the  orbitals 
(k  «  1,  2,  •••)  ,  the  above  three  functions  have  the  orthogonality 
properties  given  by  eqs.  (13a),  (13b)  and  (13c). 
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APPENDIX  II 


Equation  for  the  Difference  Function 


f(2) 


In  this  appendix  the  equation  which  the  second  order  difference 
function  f^^^  satisfies  is  developed  by  straight  forward  perturbation 
techniques. 

(2) 

The  equation  for  f  is  obtained  from  the  two-particle 
SchrCdlnger  equation 

Ht(12)  -  E+(12)  (11,1) 


by  decomposing  the  Hamiltonian  according  to  eqs.  (5)  and  (6),  expanding 
E  according  to  eq.  (1,7)  and  writing  for  \(r(12)  the  perturbation 
expansion  given  by  eq.  (12).  Equating  coefficients  of  A  and  A 
obtains,  respectively. 


where 


-  2£^j  S^*^^(12)  -  (e^^^  -  y^(l)^^(2) 


(11,2) 


(11,3) 


and 

(o) 


(■ 


-  26^1  S^^^12) 


where 


|^f^^^(l)j^(2)  +^^(l)f^^\2)J  +  -  26^j 

-  E^^^j)^(l)J?(2)  +  S^^^(12) 

(12)  H^^^^^(l)y^(2)dt^dc2  . 


(11,4) 


(11,5) 
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Multiplying  eq.  (11,4)  by  ^^(2)  ,  Integrating  over  electron  2  and 

taking  eqs.  (13a),  (13b)  aiul  (13c)  into  account,  one  obtains  the  equation 

(2) 

for  the  difference  function  f  ,  namely 

^h°(l)  -  f^^^(l)  -  ■  /s^^^*(12)H^^^^^(2)d^2  (11,6) 

where  is  given  by  eq.  (11,5)  and  satisfies  eq.  (11,2). 


